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Abstract 

The Newtonian solid-mechanical theory of non-compressional spheroidal and tor- 
sional nodeless elastic vibrations in the homogenous crust model of a quaking neu- 
tron star is developed and applied to the modal classification of the quasi-periodic 
oscillations (QPOs) of X-ray luminosity in the aftermath of giant flares in SGR 1806- 
20 and SGR 1900-1-14. A brief outline is given of Rayleigh's energy method which 
is particular efficient when computing the frequency of nodeless elastic spheroidal 
and torsional shear modes as a function of multipole degree of nodeless vibrations 
and two input parameters - the natural frequency unit of shear vibrations carrying 
information about equation of state and fractional depth of peripheral seismogenic 
layer. In so doing we discover that the dipole overtones of both spheroidal and 
torsional nodeless vibrations possess the properties of Goldstone soft modes. It is 
shown that obtained spectral formulae reproduce the early suggested identification 
of the low- frequency QPOs from the range 30 < < 200 Hz as torsional nodeless vi- 
bration modes iy{ote) of multipole degree i in the interval 2 < ^ < 12. Based on this 
identification, which is used to fix the above mentioned input parameters of derived 
spectral formulae, we compute the frequency spectrum of nodeless spheroidal elastic 
vibrations i'{qSi). Particular attention is given to the low-frequency QPOs in the 
data for SGR 1806-20 whose physical origin has been called into question. Our cal- 
culations suggest that these unspecified QPOs are due to nodeless dipole torsional 
and dipole spheroidal elastic shear vibrations: z^(oii) = 18 Hz and i^(osi) = 26 Hz. 
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1 Introduction 

The discovery of QPOs of X-ray luminosity in the aftermath of giant flares SGR 1806-20 
and SGR 1900+14 (Israel et al 2005; Strohmayer & Watts 2006; Israel 2007; Watts & 
Strohmayer 2007) with concomitant interpretation of QPOs as caused by quake-induced 
differentially rotational seismic vibrations has stimulated remarkable developments in the 
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magnetar asteroseismology (e.g. Piro 2005, Samuelsson & Andcrsson 2007a, 2007b; Lee 
2007; Levin 2007; Watts & Reddy 2007; Sotani et al 2007; Bastmkov et al 2007a and 
references therein). Following the above interpretation and presuming the dominant role 
of the elastic restoring force, the focus of most theoretical works is on computing the 
frequency spectra of odd-parity torsional mode of shear vibrations and less attention is 
paid to the even parity spheroidal elastic mode. However, from the viewpoint of modern 
global seismology (Lay & Wallace 1995; Aki & Richards 2003), the spheroidal vibrational 
mode in a solid star and planet has the same physical significance as the toroidal one 
in the sense that these two fundamental modes owe their existence to one and the same 
restoring force (e.g. McDermott, Van Horn, Hansen 1988; Bastrukov, Weber, Podgainy 
1999, Bastrukov et al 2007b). In this fight there is a possibility that, by not consid- 
ering both these modes on an equal footing, we may miss discovering certain essential 
novelties which are consequences of solid mechanical laws governing seismic vibrations of 
superdense matter of neutron stars. Adhering to this attitude and continuing the inves- 
tigations recently reported by Bastrukov et al (2007a), we derive here spectral equations 
for the frequency of both spheroidal and torsional elastic nodeless vibrations in the solid 
crust of quaking neutron star and examine what conclusions can be drawn regarding 
low-frequency QPOs whose physical nature still remain unclear. In Sec.2 by use of the 
energy variational method we derive spectral formulae for the frequencies of the nodeless 
spheroidal and torsional elastic vibrations locked in the finite-depth seismogenic layer. 
Particular attention is given to the dipole spheroidal and torsional vibrations possessing 
properties of Goldstone's soft modes. In sec.3, the obtained spectral formulae are appfied 
to a modal analysis of available data on the above mentioned QPOs. The obtained results 
are highlighted in Sec.4. 



2 Frequency of nodeless spheroidal and torsional elas- 
tic shear vibrations in homogeneous crust 

In this paper we follow the line of argument of the standard two-component, core-crust, 
model of quaking neutron star (Link, Franco, Epstein 2000) in which metal-like material 
(composed of nuclei dispersed in the sea of relativistic electrons) of the finite-depth seis- 
mogenic crust is treated as a highly robust to compressional distortions elastic continuous 
medium of a uniform density p and characterized by constant value of shear modulus 
Also, the model presumes that the quake- induced seismic vibrations driven by bulk 
force of pure shear elastic deformations, which are not accompanied by fiuctuations in 
density 5p — —p^k^k — 0, can be modeled equation of Newtonian, non-relativistic, solid 
mechanics 

PUi = VjkfTifc, (Tik = 2lIUik, Uik = \\SiUk + Vfc-Ui], Ukk = Vfe-Ufc = 0. (1) 
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From now on Ui{r,t) stands for the field of material displacements in the crust of the 
depth AR = R — R^. with R and Rc being radii of star and core, respectively. The linear 
relation between tensors of shear elastic stresses and shear deformations or strains Uik 
is the Hooke's law of elastic (reversal) shear deformations, whose strength is characterized 
by shear modulus 

In what follows we focus on poorly investigated regime of quasistatic shear vibrations 
whose mathematical analysis is quite different form vibrations in the regime of standing 
waves (Bastrukov at al 2007b). In the latter case, the solution of eigenfrequency problem 
consists in searching for the wave numbers k of standing waves of material displacements 
u obeying the vector Helmholtz equation V^u + fc^u = supplemented by two boundary 
conditions on the edges of seismogenic zone, that is, on the core-crust interface and the 
star surfacjll. The two fundamental (orthogonal and different in parity) solutions to the 
vector Helmholtz equation one of which is given by the positive parity poloidal (polar) 
field and another by negative parity toroidal (axial) field provide a basis for generally 
accepted Lamb's classification of vibrational modes as spheroidal or s-mode (in which the 
field of displacement is described by poloidal field) and torsional or t-mode (in which the 
field of displacement is described by poloidal field). The substitution of these fundamental 
solutions in boundary conditions, which are motivated by physical arguments, leads to 
the coupled system of transcendent dispersion equations for the wave vector k having 
the form of bilinear combinations of spherical Bessel and Neumann functions with highly 
involved distribution of nodes along the radial coordinate of the layer of thickness Ai?. 
It is the roots of these dispersion equations yield the discrete set of the wave numbers k 
uniquely related to the frequency of oscillations: u = kcf. The systematic application of 
the above method to the case of torsional elastic modes trapped in the homogeneous crust 
is presented in (Bastrukov et al 2007a) and pointed out that similar procedure holds true 
for the case of spheroidal vibrational mode in the homogeneous crust model. 

The situation is quite different in the case of long wavelength vibrations, that is, when 
wave vector k = {2tt/X) —>■ and, hence, the wavelength A —>■ oo. This limit corresponds 
to regime of quasistatic, substantially nodeless, vibrations in which the fields of oscillating 
material displacement subject to the vector Laplace equation (Bastrukov et al 2007b) 



Understandably, this last equation can be thought of as the long wavelength limit of 
vector Helmholtz equation describing standing-wave regime of vibrations (Bastrukov et al 
2007a; 2007b). Our interest to the regime of nodeless seismic vibrations is motivated by 

^It is appropriate to note that the very problem of vibrational modes and its solutions for the case of 
standing wave in entire volume of homogeneous elastically deformable solid sphere has first been tackled 
and solved by Lamb [11] . The application of Lamb's theory to analysis of seismic vibrations of the solid 
Earth model is extensively discussed in monographs (Jeffreys 1976; Lapwood & Usami 1981). 
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arguments of works (Samuelsson & Andersson 2007a; 2002b) in which based on equations 
of general relativity it was found that low-frequency QPOs in the above mentioned SGR's 
flare can be identified with \ow-i nodeless torsional elastic vibrations locked in the crust. 
And it is somewhat surprising that no such problem has hitherto been properly analyzed 
on the basis of equations of Newtonian, non-relativistic, solid mechanics. To this end, in 
(Bastrukov et al 2007b) it has been shown for the first time that the problem of computing 
frequency spectra of both spheroidal and torsional modes of global nodeless vibrations of 
the solid star can be uniquely solved with aid of the energy method which is due to 
Rayleigh. In the present paper this method is extended to spheroidal mode of nodeless 
elastic vibrations which are considered in one line with torsional mode. 

The stating point of the energy variational method is the integral equation of the 
energy balance 

^ / j ^ikUikdV = -2 J iiUikiiikdV (3) 

which is obtained by scalar multiplication of equation of solid mechanics, H], with Ui and 
integration over the volume of seismogenic layer. For our further purpose the field u(r, t) 
can be conveniently represented in the following separable form 

u(r,t) = a(r)a(t) (4) 

where a(r) is the field of instantaneous (time-independent) displacements obeying, as 
follows from ([2]), to equations 

V^a(r) = 0, V ■ a(r) = (5) 

and a{t) stands for the temporal amplitude of vibrations. Inserting ([3]) in ([2]) we arrive 
at equation for a{t) having the form of standard equation of linear oscillations 

dE ^ Mc? Ka^ . ■? K 

^ = 0, B = — + — - « + ^V = 0, = (6) 

M = jpa,a,dV, A- = 2/,a..a..dV a.. = i|V.a. + V.a.|. (7) 

The solenoidal fields of instantaneous material displacements in two fundamental modes 
of nodeless vibrations - the spheroidal (normally abbreviated as oSe) and the toroidal 
(abbreviated as o^^), are determined by two fundamental (orthogonal and different in 
parity) solutions to the vector Laplace equation of the vector Laplace equation built on the 
general solution to the scalar Laplace equation V^x(i") = 0- spherical coordinates with 
fixed polar axis, the solution of ([5]) corresponding to nodeless vibrations in the spheroidal 
mode, a^, is given by the even parity poloidal (polar) vector field and instantaneous 
displacements in the torsional mode, SLt, are by the odd parity toroidal (axial) vector 
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field: 



a, = V X V X (rx), at = V x (rx), (8) 
X(r) = /,(r)P,(cos^), /,(r) = [A,r' + B,r~^'^\ (9) 

Henceforth PiicosO) stands for tfie Legendre polynomial of multipole degree £ and At and 
Bi are the arbitrary constants to be eliminated from boundary conditions on the core- 
crust interface and on the star surface. Thus, the problem of computing the frequency 
spectra of both spheroidal and torsional nodeless vibrations is to fix the these constants 
and compute integrals for integral parameters of vibrations, that is, the inertia M and 
the stiffness K. 



2.1 Spheroidal mode 

The poloidal field of nodeless instantaneous displacement a^ in s-mode is irrotational: 
V X a^ = (Bastrukov et al 2007b). To specify and B(, we adopt on the core-crust 
interface, r = i?c, the condition of impenetrability of seismic perturbation in the core. 
On the star edge, r = R, we impose the condition that the radial velocity of material 
displacements equals the rate of spheroidal distortions of the star surface 



Ur\r=Ra — 0, 



Ur 



= R(t), R{t) = R[l + a{t)Pe{cose)]. 



(10) 



The solution of resultant algebraic equations leads to following values of arbitrary con- 
stants 



Ai 



i{i + iy 



B, 
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R. 
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Tedious but simple calculation of integrals for inertia M and stiffness given by ([7]), 
with poloidal field a^ yields 



M,(£, A) 



AttR^P 
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(12) 
(13) 
(14) 



The fractional frequency of nodeless spheroidal irrotational shear vibrations as a function 
of multipole degree £ is given by 



u;2(^) ^'M 2(2£+l) 



[1 - A2^+i) 



{f - 1)(1 - A2^-i) + £{£ + 2)A2^-i(l - A2^+3) 



R' 
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1^0 
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{£+l)+ ^2^+1 
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(,15) 
(16) 
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It worth noting that in the hmit of zero-size radius of the core, A = (Rc/R) 0, when 
entire volume of the star sets in vibrations, we regain the early obtained spectral formula 
for global nodeless spheroidal nodeless shear vibrations i^(os^) = z/q [2{2i +!)(£— 1)]^^^ 
showing that the lowest overtone of the global nodeless spheroidal oscillations in the entire 
volume of the star is of quadrupole degree, i = 2 (Bastrukov et al 2002a 2002b; 2007b). 
In the meantime, the lowest overtone of spheroidal vibrations trapped in the crust is 
of the dipole degree, £ = 1. This suggests that the dipole overtone can be considered 
as a signature of spheroidal vibrations locked in the crust. In the upper panel of Fig.l 
we plot the fractional frequency Us{i)/u!o as a function oi h = AR/R. This picture 
indicates that dipole vibration can be thought of as, so called, Goldstone's soft mode 
whose most conspicuous feature is that the frequency as a function of intrinsic parameter A 
of oscillating system uJs{i = 1, A) -^0, when A — > 0. In the model under consideration this 
parameter is given by A = {Rc/R)- The limit A = belongs to translation displacement 
of the center-of-mass of the star, not a vibration; this is clearly seen from the equation for 
energy (Hamiltonian) of harmonic oscillations ([6]). In the next section we show how the 
input parameters of obtained spectral equation f|T5l) . namely, the natural unit of frequency 
and the depth h of seismogenic layer can be extracted from the data on QPOs for SGS. 
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2.2 Torsion mode 



For the torsional oscillations locked in the crust, the constants Ai and Bg are eliminated 
from the following boundary conditions 



U^\r=Rc = 0, U^\r=R = [(f>R X R]<^, 



0, 



d 



d 



89 sin 9 < 



(17) 
fl8) 



First is the no-slip condition on the core-crust interface, 



Rc, implying that the 



amplitude of differentially rotational oscillations is gradually decreasing from the surface 
to the core. The boundary condition on the star surface, r = i?, is dictated by symmetry 
of the general toroidal solution of the vector Laplace equation. The support of this last 
boundary condition lends further considerations showing that it leads to correct expression 
for the moment of inertia of a rigidly rotating star. The resultant algebraic equations 
steaming from above boundary conditions lead to 



Tedious calculation of integrals for Mj and Kt leads to 

47re{e + 1) pR^ 



(2£ + l)(2£ + 3) (l-A2^+i)2 
(2£ + p2 
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{^ + 2) -2^+1 
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(21) 



:22i 



In the limit of zero-size radius of the core, A = [Rc/R) 0, corresponding to torsional 
oscillations in the entire volume of the star we regain the early obtained spectral formula 
for the global nodeless torsional elastic vibrations u^otg) = uq [(2£ + 3)(2£ — 1)]^^^ showing 
that in case of global torsional oscillations the lowest overtone is of quadrupole degree 
(Bastrukov et al 2002a; 2002b; 2007a; 2007b). However, this is not the case when we 
consider torsional nodeless oscillations locked in the seismogening layer of finite depth 
Ai? = R — Rc. For i = 1, equations (|2Ui) and ([21]) are reduced to 



Mt{i = l,X) 
Kt{e, A) = 87TfiR^ 



87rpR^ 



15(1- A3)2 
A3 



[1 - SA^ + 9A^ - 5A^] 



A3)^ 
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< A < 1. 



(23) 
(24) 
(25) 
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Figure 2: Theoretical curves for the frequency of spheroidal (dashed) and torsional 
(solid) nodeless elastic oscillations computed with aid of spectral formulae for frequency of 
spheroidal, eq.(15), and torsional, eqs. (26)-(27), modes as functions of multipole degree 
in juxtaposition with data (symbols) on QPOs for SGR 1900+14 and for SGR 1806- 
20. The modal identification is taken from (Samuelsson and Andersson 2007b; Watts &; 
Strohmayer 2007). 



In the limit zeroth core, A — > 0, the stiffness Kt ^ and the mass parameter getting 
the form of moment of inertia of absolutely rigid solid star of mass Ai and radius R: 
Mt{i = 1, A = 0) = {2/5)AiR'^. This consideration again shows that the dipole vibration 
exhibits features of the Goldstone soft mode owing its emergence to the trapping of 
torsional shear oscillations in the peripheral crust of finite depth. 

The general spectral equation for the fractional frequency of nodeless torsional oscil- 
lations of arbitrary multipole degree £, computed with aid of equations ( l20i) and ( 12T1) . can 
be presented in the following analytic form 



Pt{vo, A) = 4 



[(£ + 2)(£-l)]Mz/o,A) 



1 + 



2(£ + 2) 



'1-A 



(26) 
(27) 



1 - 



£ - A2^+i[(£ + 2) + (2£ - 1)(2£ + 3) - (2£ + lf\^ + (2£ + 



{2£ - 1) - A2^+i[(2£ - 1){2£ + 3) - {2£ + l^X^ + {2£ + 3)X^^+^] 
The practical usefulness of such representation is extensively discussed in (Bastrukov et 
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Figure 3: Theoretical predictions for the frequency of spheroidal (dashed) and torsional 
(solid) nodeless elastic oscillations as a function of multipole degree i in juxtaposition with 
data (symbols) on QPOs SGR 1806-20. The identification of QPOs pictured by squares 
is taken from (Samuelsson & Andersson 2007b; Watts & Strohmayer 2007). Based on the 
results of this latter work, our calculations suggest that low frequency QPOs discovered 
in (Israel et al 2005), that are pictured by triangles, can be identified as dipole toroidal 
and dipole spheroidal nodeless vibration, respectively: z/(oti) = 18 Hz and z/(oSi) = 26 
Hz. 

al 2007a). The fractional frequency as a function oi h = AR/R is pictured in down panel 
of Fig.l which show that the lowest overtone of torsional vibrations trapped in the crust 
is of dipole degree and that dipole overtone of differentially rotational vibrations of the 
crust against core possesses properties of the Goldstone's soft mode. 



3 Application to SGR 1900+14 and SGR 1806-20 

The obtained spectral formulae (15) and (26)-(27) describe the frequencies of both spheroidal 
and torsional nodeless oscillations as functions of the multipole degree i. The natural unit 
of frequency uq of shear elastic vibrations and the fractional depth h of peripheral seismo- 
genic layer are input parameters carrying information about material properties of neutron 
star matter (density, shear modulus) and geometrical sizes of star and seismoactive zone. 
Considering the observation data for SGR 1900+14 and SGR 1806-20, we demonstrate 
here how the obtained spectral equations can be used to eliminate some uncertainties in 
identification of QPOs. 

First, we examine the agreement of obtained spectral formula (l26|) - fl27|) for torsion 
mode with identification of the QPOs frequencies from interval 30 < u < 200 Hz with 
frequencies of nodeless torsional vibrations of multipole degree i from interval 2 < £ < 12 
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suggested in works (Samuelsson & Andersson 2007a; 2007b). In so doing we use the 
proposed in these latter works identification of QPOs in SGR 1900+14 data [namely, 
^^(0^2) = 28 Hz; 2/(0^4) = 53; Hz 1^(0^6) = 84 Hz, z/(otii) = 155 Hz borrowed from Table 1 
of paper (Samuelsson & Andersson 2007b)] as reference points and vary parameters uq and 
h entering our spectral formula (15) for torsional mode so as to attain the best fit of these 
points. The result of this procedure is shown in upper panel of Fig. 2 by solid line. Then, 
making use of the fixed in the above manner parameters uq and h, we compute (with the 
aid of spectral formula (fTSl) ) the frequency of spheroidal mode I'^oSe). The application of 
this procedure to modal analysis of QPOs data for SGR 1806-20 is pictured in down panel 
of Fig.2. Based on proposed in the above mentioned paper identification of the following 
points 2/(0^2) = 30 Hz; //(o^e) = 92 Hz and z/(otio) = 150 Hz we extract parameters z/q and 
h entering in our spectral formulae for torsional mode, equations (!26l)-(!27l). 

In Fig. 3 we highlight by triangles two non-identified before low-frequency QPOs in data 
for SGR 1806-20 (Israel et al 2005), namely z/ = 18 Hz and u = 26 Hz. The extrapolation 
to low-frequency region of spectral formulae (13) and (15) leads us to conclude that the 
latter low-frequency points can be interpreted as manifestation of the dipole toroidal, o^i, 
and the dipole spheroidal, o^i, overtones of nodeless shear elastic vibrations, respectively. 

As to high-frequency points u = 626 Hz and u = 1840 Hz in data for SGR 1806- 
20 is concerned, the lack of observational data makes the above scheme of identification 
less effective. Putting these points on spectral curve for torsional mode //(ot^) extended 
to very high value of i, we get 1^(0^^=42) = 625 Hz and z/(ot£=i22) = 1840. However, if 
one puts these points on the spectral curve for spheroidal mode i^(os^), we obtain the 
following identification z/(oS£=3o) = 625 Hz and z/(ot£=87) = 1840. The last case may be 
more favorable because the higher multipole degree of oscillations the less their lifetime 
(Bastrukov et al 2007b), and low-£ overtones have, therefore, more chances for surviving. 
But one must admit that this viewpoint is highly questionable and should be thought of 
as suggestive, not conclusive. 

4 Summary 

The exact spectral formulae, which has been obtained here within the framework of New- 
tonian, non-relativistic, solid-mechanical theory of seismic vibration for the first time, are 
interesting in its own right from the viewpoint of general theoretical seismology (e.g. Lay, 
Wallace 1995) because they can be utilized in the study of seismic vibrations of more wider 
class of solid celestial objects such as Earth-like planets. One of the remarkable findings of 
our investigation is that the dipole overtones of nodeless elastic shear vibrations trapped 
in the finite-depth crust of seismically active neutron star possess properties of Goldstone 
soft modes. It is shown that obtained spectral equations are consistent with the existence 
treatment of low-frequency QPOs in the X-ray luminosity of flares SGR 1900+14 and 
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SGR 1806-20 as caused by quake-induced torsional nodeless vibrations (Samuelsson & 
Andersson 2007a; 2007b). What is newly disclosed here is that previously non-identified 
low-frequency QPOs in data for SGR 1806-20 can be attributed to nodeless dipole tor- 
sional and spheroidal vibrations, namely, z^(o^i) = 18 Hz and i^(oSi) = 26 Hz. 
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